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Euclidean point-to-point propagators or wall-to-wall correlators related to exchange by an unstable 
particle (a, p, tj-mesons) are modified by presence of particle width. In particular, the usual method 
of deriving particle masses from logarithmic derivatives need to be modified. Similarly Yukawa- 
like potentials of nuclear physics due to exchange of those mesons are significantly modified since 
the coupling to the decay products is strong. For example, the large distance asymptotic changes, 
~ exp(— M m inT), where M m i n is the sum of the decay product masses (2m- n , 2m e , 2m„). In the area 
M m i n r < 1 the potential has a long-range tail ~ 1/r 3 . Similar effects appear due to the virtual 
decays in the elecroweak sector of the Standard model. The Z — 7 mixing via electron loop gives 
the parity violation potential with the range l/2m e , i.e. the range of the weak interaction increases 
10 s times. 

PACS numbers: 21.30-x,12.38.Gc,21.30.Cb,13.75.Cs 



I. INTRODUCTION 

This paper addresses the following general question: 
what is the role, if any, of a width T of an unstable par- 
ticle in situations in which this particle is a virtual in- 
termediate state so that its physical decay cannot take 
place. We will discuss three applications of such kind: 

(i) Euclidean point — to — point correlation functions in 
which particles in question are emitted and absorbed by 
two operators acting at 4-d points; 

(ii) Euclidean wall — to — wall correlators, widely used 
on the lattice for determination of particle masses; in this 
case a signal is emitted and absorbed from 3-d planes sep- 
arated by some time interval r; 

(iii) Yukawa-like potentials due to meson exchanges 
widely used in nuclear physics, this case can be called 
line — to — line correlators. 

The cases (ii) and (iii) are schematically illustrated in 
FigJU although those have direct applications, the corre- 
sponding results simply follow from those of the generic 
case (i), the point-to-point Euclidean propagator. There 
are of course other similar situations in which the same 
ideas apply but not discussed below: t-channel meson 
exchange in scattering amplitudes is another example. 
There would be no difficulty to generalize what is done 
below to these case as well. Also let us add that al- 
though we call unstable particle a meson, it can equally 
well be any other unstable particle with non-negligible 
width (e.g. the Z boson of weak interaction). 

To set up a problem, let us first explain why naive Eu- 
clidean continuation of the real-time amplitudes does not 
work. Usual real-time amplitudes describing resonances 
behave as damped oscillators 

A(t) ~ exp(-iMt - Tt/2) (1) 

Their naive Euclidean rotation, by simple substitution 



FIG. 1: The setting for (a) wall-to-wall and (b) line-to-like 
Euclidean correlators, describing Euclidean lattice correlators 
and exchange potentials, respectively. The unstable article 
propagation is indicated by a dashed line. 



t — ► — it which leads to an oscillating function ~ 
exp(—Mr + iTr/2) which certainly cannot be true Eu- 
clidean propagators or amplitudes. In fact those are de- 
fined via real exponent of the Hamiltonian ~ exp(-Hr) 
and standard spectral decomposition plus Hermitian 
Hamiltonian of course demands all energy to be real and 
thus the correlator to decrease monotonously with r. 

Why should one seek an answer to such questions? 
Apart of being the basis of applications of statistical me- 
chanics, Euclidean correlation functions with real expo- 
nentials of H are the main tool of lattice gauge theory 
community, which studies the non-perturbative dynam- 
ics of non-Abelian gauge fields by numerical simulations, 
see multiple reviews e.g. [l[ . Rotation to Euclidean time 
is also a crucial ingredient of the semiclassical approach 
based on instantons, see e.g.reviewQ- These correlators 
are widely used to extract the particle masses and their 
couplings to certain currents, but it was not yet used to 
extract the widths of resonances [Hj]. We will see below 
that it is indeed not a simple task. 

The present-day lattice calculations have practical lim- 
itations stemming from available computer resources. As 
a result, quark used are still rather heavy, so that in most 
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calculations the pions are so heavy that no decays of a 
and p are possible. Another kinematical effects is a conse- 
quence of small size of the present-day lattices L ~ 2 fm. 
Indeed, momentum quantization of the decay products 
due to the periodic boundary conditions to p n = 2im/L 
also prevents decay in many cases, especially with the 
the initial "plane-to-plane" projection[13| to p = 0. 

This paper is not addressing these lattice limitations, 
considering ideal case with quark masses at its physical 
value and L large enough to cause no problem. The is- 
sue we address is how a propagator of unstable particle 
look like, when it is virtual and is propagating either in 
space direction (related to Yukawa potentials below) or 
in imaginary (Euclidean) time. As we will see, the widths 
can modify the shape of these propagators , sometimes 
beyond recognition. Although masses and widths can 
still be extracted from those correlators, it would require 
very complicated fitting procedure. For the same reason, 
using lattice simulations for the evaluation of the trans- 
port properties via Kubo-like formulae from Euclidean 
correlators is practically impossible, as discussed e.g. in 

1- 

Our approach is based on standard dispersion rela- 
tions, which connect real part of the amplitudes in mo- 
mentum representation to their imaginary part, also 
known as spectral densities, in which resonances appear 
as Breit-Wigner-like peaks. If the resonance width is 
very small, those can be approximated by delta func- 
tions and the width ignored. If so, the usual exponen- 
tial behavior of the wall-to-wall Euclidean correlators 
K(t) ~ exp(—Mr) or standard Yukawa-like potentials 
V(r) ~ exp(—Mr)/r are obtained. If the width is small, 
one can perform a saddle point integration and find an- 
alytically the first corrections to delta function, in form 
of some "shifted mass" depending on the width T. If the 
width is not small, the form of the correlators or poten- 
tials is strongly modified, as it will be detailed below. 

Another related issue is modification of hadronic reso- 
nances at finite temperatures. We will not discuss here 
multiple physical effects which are important in reso- 
nance emission from "heat baths" created e.g. in heavy 
ion collisions, but just comment on the simplest "kine- 
matical" effect [U which appear even in high energy ele- 
mentary pp collisions. When the spectral density is con- 
voluted with the thermal Boltzmann exponent, the finite 
width of a resonance leads to different weighting of the 
left and right side of the peak and thus an apparent reso- 
nance mass shift. For small width and/or high tempera- 
ture it is very easy to estimate the effect: approximating 
the resonance as 

exp(-M/T) 
l + (4/r2)(Af-m res )2 ~ {) 



exp[-(-^){M 2 - 2Mm res - M 2 + 



one gets the mass shift to be 

r 2 

Sm max « mo - — (3) 

ol 

This simple analytic expression however does not work 
well in practice and real shift depends much on the res- 
onance shape (see Q for details). 

II. POINT-TO-POINT CORRELATORS, 
DISPERSION RELATIONS AND THE 
SPECTRAL DENSITY 

We start with the point-to-point correlation functions 
because other cases can easily be obtained from them, 
by integration. In fact point-to-point correlators in Eu- 
clidean space-time are not just methodical object, they 
are important tools widely used in studies of structure 
of the QCD vacuum. They can be deduced phenomeno- 
logically, using vast set of data accumulated in hadronic 
physics. Second, they can be directly calculated ab ini- 
tio using quantum field theory methods, such as lat- 
tice gauge theory, or semiclassical methods. Significant 
amount of work has also been done in order to under- 
stand their small-distance behavior, based on the Op- 
erator Product Expansion, see 0] and vast subsequent 
literature. 

These correlation functions are vacuum expectation 
value of the product of two (or more) of operators 

K(x-y)=<0\O(x)O(y)\Q> (4) 

which can have any gauge invariant combination of fields, 
for example mesonic operators we will consider have local 
quark-antiquark combination 

O mes (x) = ipiMijipj (5) 

where a matrix M can include various flavor and Dirac 
structures, and a convolution over quark color indices. 
Since the vacuum is homogeneous, the correlators depend 
on the relative distance only. If two points are inside the 
light cone, a lot of things may happen: e.g. a meson may 
decay in between. We will not consider such case below, 
always assuming that the distance between the points 
(x-y) is space-like. It may really be a space distance or 
"Euclidean time" interval: in both cases there are no 
real decays and one deals with monotonously decreasing 
functions of the distance. 

The Fourier transform of K{x), K mom (q ), depends 
on the momentum transfer q flowing from one opera- 
tor to another. In its language, the situation we con- 
sider corresponds only to virtual space-like 4-momenta q 
Q 2 = —q 2 > 0, like in scattering experiments in which 
the resonance is in t and not s channels. 

Due to causality, standard dispersion relations follow: 
they relate real and imaginary parts 

2 f dsImK mom {s) 

K mom {q = / — (6) 

J n (s - q z ) 
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where the r.h.s. is the physical spectral density 
ImK mom {s). It certainly is non-zero only for positive s 
above certain threshold (twice the pjon mass in the case 
of p, a resonances to be considered jl4{. 

Dispersion relation is the basis of various sum rules. 
Their general idea is as follows; suppose one knows 
the l.h.s. K mom {q 2 ) somewhere; thus some integral in 
the r.h.s. of the physical spectral density is known as 
well. (For example, for mesonic operators the correla- 
tors at small x are just given by free quark propagation 
K ~ with calculable corrections.) The simplest 

of them use directly momentum space: those are known 
as finite energy sum rules. Unfortunately, those are not 
very useful because most of the dispersion integrals are 
divergent and to usable sum rules appear only after some 
"subtractions" , which introduces extra parameters and 
undermine their prediction power. 

Improved sum rules appear if one takes a sufficient 
number of derivatives of the dispersion relation at Q = 0, 
defining the so called moments of the spectral density 



ds ImKr, 



(7) 



Following Shifman,Vainshtein and Zakharov this 
method is traditionally used in the discussion of "char- 
monium sum rules" related to correlators of cc currents. 
Another idea also suggested by these authors [5( is to 
introduce the so called Borel transform of the function 
K mom (Q) defined as follows 

K bor (m) = lim nQ 2^ oom 2 = Q2/ n 2 =fixed 



Q 2 



(-^o) n K mom (Q 2 ) (8) 



(n-l)P dQ 2 
leading to "Borel-transformed" sum rules 

/ds s 
— ImK mom {s)exp{ -) 
7r m z 



(9) 



in which the integral is cut off at large s exponentially. 
Such form of the sum rules have been widely used in 
multiple papers based on the QCD sum rules. 

Another natural option, advocated by one of us in Q, 
is to use the sum rules directly in the (Euclidean) coordi- 
nate representation. By applying Fourier transformation 
to © one obtains the following nearly self-explanatory 
form 



K(x) 



ds 



ImK mom (s)D(s 1/2 ,x) 



(10) 



The former function in the r.h.s. describes the amplitude 
of production of all intermediate states of mass \fs, while 
the latter 



D(m,x) — (m/4:ir 2 x)Ki(mx) 



(11) 



is nothing else but the Euclidean propagator of a mass-s 
state to the distance x. At large x the propagator goes 



as exp(-mx), thus it is also exponentially cut off and, 
in practice, there is little difference between this expres- 
sion and Borel sum rules . However the coordinate form 
can be much easier evaluated numerically, on the lattice 
or in scmiclassical models, unlike the Borel-transformed 
version. As discussed in [6j and subsequent literature, for 
vector and axial correlators one can directly use exper- 
imental data and calculate the l.h.s. with rather small 
error bars. This can provide a very good test of lattice 
calculations. The inverse transformation - from correla- 
tors to spectral densities - is a badly defined problem, 
quite difficult in practice, except in some simplest situa- 
tions. 

Now, let us think about a situation in which the lowest 
state in the corresponding channel is an unstable particle 
- for example p, lo mesons for vector operators or a for 
scalar isoscalar operator qq. In the simplest case of very 
narrow resonance in which the width of the resonance is 
neglected, one simply gets the propagator of a state with 
the resonance mass M K(x) ~ D(M, x). If the width is 
small, one can use the same idea as for the thermal case, 
leading to analog of the formula ([3]) with the temperature 
T replaced by l/x. If the width is substantial, the inte- 
gral in (fTU)) should be calculated more accurately, with 
the appropriate shape of the resonance. We will show 
how this procedure works in the next section. 



III. WALL-TO-WALL CORRELATORS AND 
LATTICE MEASUREMENTS OF HADRONIC 
MASSES 

Although there is quite substantial literature on lattice 
point-to-point correlators, most of studies rely on wall- 
to-wall correlation functions, obtained from K(x) by ad- 
ditional integration over the 3-dimensional spacial plane. 
This puts the momentum of all intermediate states to be 
zero, so dispersion relation in this case can be considered 
to be related not to the mass of the state but simply to its 
energy. The resulting function is related to the physical 
spectral density by 

/dm 
ImK mom (m)exp(-Tm) (12) 

where 4-d propagator D is substituted by its 1-d version. 

The mass of the lowest hadrons is traditionally ob- 
tained from the logarithmic derivative of such correlators 



E{t) 



dlogK 
dr 



M 



(13) 



Indeed, multiple lattice calculations have been able to 
show that the logarithmic derivative of wall-to-wall cor- 
relators do have "plateaus" independent of r, being the 
standard source of the meson masses. 

Such procedure works well in present-day calculations 
in which kinematical effects mentioned above prevent de- 
cays. However, as the algorithms and supercomputers 
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FIG. 2: (a) The shape of p (right) and a resonance, as a func- 
tion of its mass m (/m _1 ). (b) The corresponding logarithmic 
derivative of the wall-to-wall correlator E(t) (GeV) vs Eu- 
clidean time r(/m). The solid curves are for p (higher) and a 
resonances (both mass and width are 450, MeV), two dashed 
curves are for p with the width reduced to To = 15, 1.5 MeV, 
(lower and upper, respectively). 



gets more powerful, the widths eventually would become 
an issue. 

Let us demonstrate what should happen with "mass 
plateaus" by calculating them for different resonance 
shape. The simplest case in which an analytic formula is 
easily obtained is a Gaussian resonance ~ exp[— 4* (m — 
77i0) 2 /r 2 ] for which the mass integral can be taken from 
minus to plus infinity and 



E{t) 



m — tT 



(14) 



But this dependence does not reproduce what happens 
for realistic shapes of the resonances, with a particular 
dependence of the width on energy, especially near the 
threshold. The best way is to calculate numerically the 



behavior of E[t) for realistic spectral densities, and see 
what should its true shape for vector and scalar currents 
be. The most studied vector currents have two lowest 
resonances, 1=1 p meson and 1=0 u>. In principle, the 
spectral densities describe continuum of states, all the 
way from the lowest physical states (2 and 3 pions at rest) 
respectively. So, one may only argue that the logarithmic 
derivatives is limited by 2771^ or Sm^ from below. 

A reasonably good parameterization of the spectral 
density ImK mom for binary channels is Brcit-Wigner 



p(M) = 



r(M) 



(M - m res y + T(M) 2 /4 



(15) 



where the width is not constant and depends on the cur- 
rent running mass M and angular momentum of the re- 
sulting pions. Due to P-wave nature of p — > 7T7t decay, 
its r p (M) = 4rn 2 /M 2 ) 3 / 2 , while for s-wave a res- 

onance one should use the power 1/2 instead. The cor- 
responding shape of the resonances is shown in FigJ^a): 
we do not show here fitting of the real data, although 
we have checked that these shapes do indeed agree with 
them. 

The resulting "realistic" logarithmic derivative of the 
meson contribution is shown in Fig[2]Jb). Let us start its 
discussion with two dashed curves, which are for p me- 
son with the width artificially reduced by the factor 100 
(upper dashed) and 10 (lower dashed curve) compared to 
real- world To = .14 GeV. These two curves nicely show 
"plateaus" at the rho mass M p = .77 GeV at not-too- 
large r, with a "slide" toward the threshold (Im^) at 
large r. However for the realistic width value (the solid 
curve) the rho meson mass plateau is already gone, with 
a slide of E(t) even at small t. The situation is even 
more dramatic for the sigma meson, for which the mass- 
to- width ratio is 1:1. Presumably by fitting such curves, 
if they come from correlator studues, one would still be 
able to recover both the mass and the meson width, but 
we leave this for further studies. 

It may be useful to have some analytical estimates for 
K (t). Let us start from large r. In this case small m 
dominate in the integral (|12[) and we can substitute m = 
2m w (l + x) assuming x « 1. The result is 



K(t) = const 



E(t) = 



2m 1l To exp{—2m^r) 
(m res - 2m 7r ) 2 (m w r)« +1 
dlogK 



dr 



2m, + (g + l)/r, 



(16) 
(17) 



where q = 3/2 for p and q — 1/2 for a. This should 
be compared to the resonance contribution K r {r) — 
exp(—m r T). Only if this term dominates, we have 
E t (t) « m r . However, for p and a the coefficient 
2m 7V To/ (m res — 2m 7r ) 2 is not small, therefore, the asymp- 
totic regime starts very early. The transitional area corre- 
sponds to 2m- K < 1/t < m r . Here an essential contribu- 
tion comes from the integration interval 2m n < m < m r : 
K tr (T) ~ Tq/ '(ir(m res — 2m T ) 2 r). If this contribution 
dominates we have E tr {r) w 1/r. Note that the value of 
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the width actually determines the position of the transi- 
tion point between the areas of the resonance and non- 
resonance domination in K{t) or E(r), Therefore, both 
the meson mass and width can, in principle, be found 
from the fit of lattice calculations. 



IV. YUKAWA-LIKE POTENTIAL DUE TO 
EXCHANGE BY AN UNSTABLE PARTICLE 

Let us consider the potential between two nucleons pro- 
duced by a meson (e.g. a, p, u>) exchange. For brevity we 
omit all isotopic and kinematic structures which are dif- 
ferent for different particles and may be easily restored in 
the final answers. This potential in the coordinate repre- 
sentation may be presented as an integral of its imaginary 
part in the momentum representation 



U(r) 



1 



47r 2 r 



ImU(M) exp(-rM)2MdM , (18) 



One can either consider this expression as simply a line- 
to-line variant of the Euclidean correlators, with 3-d 
Green function exp(—Mr)/r substituting 4-d one in (TIT))) , 
or follow the usual derivation of the Uchling potential in 
textbooks (such as 0). The lower limit of the integra- 
tion is determined by the minimal mass 2m of the decay 
products (7r-mesons) of the intermediate particle and 



ImU(M) 



Tm 



4tt5i3 2 



M 2 - M 2 - P{M) 
A-Kg 1 g 2 m r T{M) 
(M 2 - Ml) 2 + m 2 T 2 (M) 



(19) 
(20) 



Here (71,(72 are the meson- nucleon interaction constants, 
P(M) is the polarization operator, Mo is the bare meson 
mass, M m = M m (M) = 

= Mr, 



v/M 2 + ReP(M) is the running 



meson mass, m r — M m (m r ) = yj M$ + ReP(m r ) is the 
meson mass at the resonance, and ImP(M) = — m r T{M) 
where T(M) is the running meson width. 

Substituting Eq. JTSJ into Eq.(HHJ) we obtain 



U(r) 



g\g%m r 



T(M) exp(-rM) 
(M 2 -M^) 2 + m 2 r 2 (M) 



dM 2 (2\) 



The explanation of this formula, like other Euclidean 
correlators above, is simple: the mass M of the unstable 
particle is not fixed and therefore the potential can be 
seen as an the integral over all possible meson masses 
given by the spectral density (e.g. a modified Brcit- 
Wigner formula we use). 

In the case of vanishingly small T(M) the Breit-Wigner 
factor in Eq. (|21[) is proportional to ^-function and the 
integration gives the usual Yukawa potential: 



U (r) = cxp(— rm r ) 



(22) 



Non-zero width T(M) produces significant deviations 
from the simple Yukawa formula. The potential U(r) in 



4 

3.5 

* 3 

3 

03 

M 2.5 
X 

i 2 
O 1.5 
1 

0.5 





/ 

/ / 




/ / 




// 




// 



















FIG. 3: The ratios of p-meson potential to Yukawa potential 
(thick solid line) and a-meson potential to Yukawa potential 
(dashed line) as functions of m r r. At m r r < 4 the negative 
corrections given in Eq, (130[l dominate. At m r r > 4 the pos- 
itive corrections Eqs. (| 2411 26\ change the long-range behavior 
of the potentials. 



Eq. (|2ip can be easily found using numerical integration. 
At Fig[3]we presented graphs of two ratios: p-meson po- 
tential divided by the Yukawa potential and er-meson po- 
tential divided by the Yukawa potential . We see that the 
deviations from the Yukawa formula may be very large, 
especially at large distances. Below we perform analyti- 
cal estimates to explain the nature of different corrections 
to Yukawa formula. For the analytical estimates we as- 
sume m << m r and T r << m r where T r = T(m r ) is the 
resonance width. 

The integral has four important intervals of the inte- 
gration: the first interval is near the threshold, M ~ 2m; 
the second interval is 2m << M << m r ; the third in- 
terval is near the resonance, M ~ m r ; and the fourth 
interval is T r << M — m r < infinity. Therefore, the 
potential may be presented as a sum of four terms 



U(r) = U 2m (r) + U p (r) + U res (r) + U inf (r) 



(23) 



Consider first a behavior of the potential (|2lj) at large dis- 
tance, rrn r >> 1. In this case the threshold interval M ~ 
2m gives the dominating contribution. The width r(M) 
vanishes near the threshold as T(M) = (1— 4m 2 /M 2 ) q -T 
where q — 1/2 for cr-meson (the 7r-mesons are in s-wave; 
for Z-boson q = 1/2 too) and q = 3/2 for p- meson (the 
7r-mesons are in p-wave). Taking M = 2m (1 + x) << m r 
we obtain 



U2m{r) 



D q = 



D q gi92 exp(-2mr) 
r(mr)i +1 

m 2 r n 



2 3 / 2 -i(2 q y.\m r m 2 T 
(M^(2m) — Am 2 ) 2 m r (m r — 2m) 5 



(24) 
(25) 



In the estimate we assumed M m (2m) w m r . This contri- 
bution has the same nature as the Uehling potential pro- 
duced by virtual "decay" of photon to electron-positron 
pair (q = 1/2). 

Consider now the interval of the integration 2m << 
M << m r which corresponds to the distances l/m r << 
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r << l/m. The result of the integration in this area is 

2g 1 g 2 T{r- 1 ) 



U P {r) 



(26) 



The potentials U p oc 1/r 3 (l/m r << r << l/m) and 
U% m oc exp(— 2mr)/r q+2 (r >> l/m) cover the area 
r >> l/m r . The potential here is exponentially en- 
hanced in comparison with the Yukawa potential- see 
FiglU In the region r ~ 1/ m r the resonance contribution 
dominates. To estimate the correction to the Yukawa po- 
tential we present the integral in Eq. (|2Tj) in the following 
form: 



= g l92 m r eM-rm r ) I{r) 



I(r) = 



exp[r(m r - M)]T r 2MdM 
, m [(M 2 -to 2 ) 2 +to 2 F 2 ] 
[1 - r{M - m r )}T r 2MdM 



(27) 
(28) 
(29) 



/n [(M 2 - to 2 ) 2 + m 2 T 2 } 

This gives an estimate of the corrections to the Yukawa 
potential 

9i 92 



U res {r) 



■ exp(— rm r )(l — T r r/ir — T r /(nm r )).(30) 



Thus, to use Yukawa potential we must assume T r << 
Trm r and T r r << ir. Then we should add U p or U2 m m 
the area m r r > 1. 

The behaviour of the potential at small distances de- 
pends on the convergence of the integral in the area 
to << M < 1/r. For mesons of finite size the integral 
is convergent at r = and Ui n f{r) should be neglected. 
The meson exchange theory of strong interaction is not 
applicable at small distance anyway. 

Now we can compare our analytical estimates with nu- 
merical calculations presented at FigJ3] To calculate the 
integral in Eq.(f2T|) we take 
r(M) = [(1 - 4m 2 /M 2 )/(l - 4to 2 /to 2 )]9 • r r ; 
T r /m r = 0.19, 2m/m r — 0.36 and q — 3/2 for p-meson; 
r r /m r = 1, 2m/m r — 0.5 and q = 1/2 for cr-meson. The 
condition T r << m r is not satisfied for cr-meson. Nev- 
ertheless, the negative corrections 1 — T r r/ir — T r / '(jrm r ) 
given in Eq. (|30p correctly explain behavior of both po- 
tentials at m r r < 1 (including the negative slope for 
m r r < 0.5, see FigU). For m r r < 3 the p-meson po- 
tential in Eq. pTj) is about 10% smaller than the Yukawa 
potenatial, the tr-meson potential is about 40% smaller. 
At m r r > 5 the positive corrections Eqs. (|24H2"6")) give 
dominating contribution and change the long-range be- 
havior of the potentials [15j. 

Finally, let us consider potential U(r) in the short- 
range approximation, 

,4^51^2 



U(r) « B- 



-S(r), 



B = 



^9\9i 



U(r)d 3 r 



r ){x)dx 



2m/m, 



x[(x 2 



(31) 
(32) 
(33) 



\ 



2 
m r r 



FIG. 4: The ratios of p-meson potential to Yukawa potential 
(thick solid line) and cr-meson potential to Yukawa potential 
(thick dashed line) as functions of m r r. Thin dashed lines 
show factors 1 — T r / (nm r ). 



where x = M/m r and "f(x) = T(M)/m r . This short- 
range approximation may also be useful in the calculation 
of the nuclear mean field. The numerical integration gives 
B = 0.93 for p-meson and B = 0.62 for cr-meson. 

It may also be useful to have some analytical estimates 
for B. In the logarithmic approximation 



B 



r , to 

(In 



Am 2 



st)]. 



(34) 



It is interesting that the logarithmic correction diverges 
at to = 0. This happens because of the long-range char- 
acter of the potential in the area l/m r << r « l/m. 
Indeed, the integration of U p (r) oc 1/r 3 from Eq. (j2"6|) im- 
mediately gives this logarithmic correction to B . How- 
ever, for real mesons ln[(m r /(2m)] ~ 1 and the negative 
correction factor 1 — T r /(Trm r ) to the resonance contribu- 
tion ( see Eg. pop ) is more important. In fact, numerical 
results for B are surprisingly close to this correction. In- 
deed, 1 — r r /(nm r ) — 0.94 and B — 0.93 for p-meson; 
1 — r r /(7rm r ) = 0.68 and B = 0.62 for cr-meson. This 
gives a very simple rough estimate for the effect of large 
width for the mesons: multiply the Yukawa potential by 
the factor 1 — T r / (■nm r ) . The results of this prescription 
in comparison with real potentials are shown at Fig[5] 



LONG-RANGE PARITY VIOLATING 
POTENTIAL 



The Z-boson exchange creates a parity violating po- 
tential. A possibility of a longer range of such poten- 
tial could have interesting applications, for example, in 
atomic tests of the Standard Model (see e.g. review of 
atomic parity violation in [9J] and a different mechanism 
for a long-range parity violating potential in [l(|). Note 
that the case of Z-boson differs from the meson case. The 
imaginary part of the polarization operator for point-like 
particles like photon and Z-boson rapidly increases with 
M. For example, the electron contribution to the photon 



7 



vacuum polarization is [7] 



ImP(M) 



-|(l-^)V 2(2m2+ M 2 ) (35) 



For the Z-boson we need to sum over all decay chan- 
nels, however, it does not change the conclusions. First 
consequence of ImP(M) cx M 2 is that the high-energy 
contribution to the integral (|21|) gives famous running 
coupling constants (correction Ui n f(r) cx In mr/r). The 
second consequence is the strong dependence of the run- 
ning width r(r _1 ) on r _1 in Eq. (f2u| since in this case 
r(r _1 ) ~ a(m 2 + r~ 2 )/m r . Therefore, the dominating 
contribution in the area l/m r << r << 1/m is actu- 
ally U p (r) oc l/M^r 5 . As a result, the correction to the 
short-range approximation, U(r) — ^V 92 5(r), remains 
finite for m = (for example, for zero neutrino mass). 
However, the Z — 7 mixing via electron loop gives a long- 
range contribution (the second term below) to the parity 
violating electron-nucleus interaction 

w G r n M j_ v 2 ^ 1 Z 4sin 2 9 w )m 2 

= 2Vl 7 ~ 3^ /( - r ^ 36 ) 

f°° 1 

I(r) = / exp(-2xm e r)y/x 2 - 1(1 + —=■)(&: (37) 
A 2a; 2 

Here the first term is the usual parity violating electron- 
nucleus interaction (see e.g. @); G is the Fermi constant, 



75 is the Dirac matrix, Q w w — N + Z(l — 4:sin 2 6 w ) is 
the weak charge, Z and N are numbers of protons and 
neutrons, 9 W is the Weinberg angle, p(r) is the nuclear 
density normalised to 1 (p(r) ~ 5(r)), m e is the electron 
mass. To obtain this potential we used the imaginary 
part of the polarization operator from Eg. ([55]) and the 
Standard model values of the interaction constants [l6l |. 
The integral has the following asymptotics: for m e r <C 1 



I(r) 



(2m e r) 2 



(38) 



This gives W cx 1/r 3 . For m e r 1 



J(r) 



/7F exp(— 2m e r) 
4(m e r) 3 / 2 



(39) 



This gives W cx exp(— 2m e r)/r 5 / 2 . Thus, we have the 
parity violating potential with the range l/(2m e ). A sim- 
ilar effect appears for the nuclear-spin-dependent part of 
the parity violating interaction. 
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